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Quasi-Concavity for Gaussian Multicast Relay 

Channels 

Mohit Thakur and Gerhard Kramer 


Abstract —Upper and lower bounds on the capacity of Gaussian 
multicast relay channels are shown to be quasi-concave in 
the receiver signal-to-noise ratios and the transmit correlation 
coefficient. The bounds considered are the cut-set bound, decode¬ 
forward (DF) rates, and quantize-forward rates. The DF rates are 
shown to be quasi-concave in the relay position and this property 
is used to optimize the relay position for several networks. 

Index Terms —capacity, decode-forward, multicast, relaying 


I. Introduction 

A multicast relay channel (MRC) is an information network 
with a source node, a relay node, and two or more destination 
nodes, and where one message originating at the source 
should be received reliably at the destinations. We consider 
additive white Gaussian noise (AWGN) MRCs and show that 
certain information rate expressions are quasi-concave in the 
receiver signal-to-noise ratios (SNRs), the transmit correlation 
coefficient, and the relay position. Quasi-concavity means 
that efficient algorithms can optimize signaling and the relay 
position. For example, suppose one wishes to place a relay to 
maximize the multicast rate ID> We studied this problem 
for a decode-forward (DF) strategy in the low-SNR regime in 
@-0. This paper extends the results to other strategies such 
as quantize-forward (QF) and to general SNR regimes. 

The paper is organized as follows. Section |I^ presents 
the MRC model and bounds on the MRC capacity. Sec. m 
shows that these bounds are quasi-concave in a correlation 
coefficient and the SNRs. Sec. HYl shows that the best DF rate 
is quasi-concave in the relay position. Sec. |V] compares the 
performance of different DF strategies, and computes optimal 


relay positions. Sec. VI discusses complex channel models and 
Sec. |VII| concludes the paper. The Appendices provide useful 
results on concavity and quasi-concavity. 


II. MRC Model and Information Rates 
A. Model 

An MRC has three types of nodes: 

• a source node s that generates a message W and transmits 
the symbols X” = X,s^i,Xs, 2 , ■ ■ ■, 
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Fig. 1: MRC with two destinations. 


• a relay node r that receives and forwards symbols Yr^k 
and Xr,k, respectively, for fc = 1, 2, • • • , n; 

• destination nodes j = 1, 2,..., W where node j receives 
Yp = Yyi, Yj^ 2 , ■ ■ ■, Yj,n and estimates W as Wj. 

We denote the destination node set as T = {1, 2,..., W}. The 
classic relay channel has = 1 and Fig. [T] shows an MRC 
with N = 2. A memoryless MRC has a function h{-) and 
a noise random variable Z so that for every time instant the 
N + 1 channel outputs Y = (Yr Yi ... Yv) are given by 

Y = hiXs,Xr,Z). 


The noise Z is statistically independent of Xs and Xr, and the 
noise variables at different times are statistically independent. 
An encoding strategy has 

• a message W uniformly distributed over {1,2,..., M}; 

• an encoding function Cg such that AT" = es(IV); 

• relay functions Cr^k with X^ k = er,/c(Y.,i,.., Y,fc-i), 
where fc = 1, 2,..., n; 

• decoding functions dj such that dj{Yp) = Wj, j G T. 


The error probability at destination j is Pe,j = Pr ^ Yl 
The multicast rate is R — (log 2 M)/n bits/use. The rate R is 
achievable if, for any e > 0 and sufficiently large n, there is 
an encoding strategy with Pej < e for all j G T. The capacity 
C is the supremum of the achievable rates. 


B. Information Rates 

The following bounds were given in |j^ for the relay channel 
{N = 1). Their extensions to MRCs are straightforward. 

• Cut-Set Bound: C < Res where 


Rc 9 = max < min 

\l<0<N 

min(/(YgY,; Y), I(X,-, YrYfXr))} (1) 
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Fig. 2; AWGN-MRC with two destinations. 


and where the maximization is over all XgXr- 

• Direct-Transmission (DT) Rate: C > Rdt where 

Rdt = max{ min liXg'^YAXr = a;*)} (2) 

l<i<Ar ■' 

and where the maximization is over all x* and Xg- 

• DF Rate: C > Rdf where 

Rdf = max < min 
y^<3<N 

xm^{I{XgXr] Yj), I(Xg-Yr\Xr))} ( 3 ) 

and where the maximization is over all XgXr- 

• QF Rate: C > Rqf where 

i?oc’=max< min min (/(XsATj.; K.) — 

/(F,; %\XgXrY^),IiXg-XY^ |X,)) } (4) 

and where the maximization is over all XgXXr such that 
Xg and Xr are independent and Xg — XrYr — Y^ forms 
a Markov chain. 

C. AWGN MRC 

The (real) AWGN-MRC has real channel symbols and 


We consider individual average block power constraints 


E 




.fc = l 


< nPg, E 




.fc=i 


< nPr- 


( 9 ) 


The SNR and the capacity of the link from node u (with 
transmit power Pu) to node v are the respective 


C(SNR„,,) = ilog(l + SNR„.,). 

We simplify the above rate bounds for the AWGN-MRC. 
• Cut-Set Bound: 


( 10 ) 

( 11 ) 


Res = max 

p 


min min (CfSNRc 1 + SNRr 1 + 

l<j<N V ' J 


2pv/SNR~SN^) , C((l - p2)(SNR«j -I- SNR«,,)) 


( 12 ) 


where p = E [XgXj.] /a/E [AT^] E [X'^] satisfies \p\ < 1. 
One can restrict attention to non-negative p. 

• DT Rate: 


• DF Rate: 


Rdt = min C(SNRs ,). 

i<j<N 


( 13 ) 


Rdf = max 

p 


min min (C(SNRs + SNRj. ,-|- 
l<j<Ar ' V .J .J 


2pXsm~sm~), c((i - pX^Rs^)) 


(14) 


One can again restrict attention to non-negative p. 

• QF Rate: Optimizing XgX^Yj. seems difficult. Instead, 
we choose Xg and Xr to be zero-mean Gaussian with 
variances Pg and Pr, respectively. We further choose 
Yr = Yr + Zr where Zr is zero-mean Gaussian with 
variance Nr- Optimizing gives (see ||^ p. 336-337]) 


Rof = min C(SNRs,-l- 
l<j<N V 


SNR^jSNR, 


SNR 




SNR, 


r,] 


SNRjj r + 1 


( 15 ) 


Yr — (Xg rXg -\- Zr ( 5 ) 

Yj — Clg jXg “t“ Clr^jXr Y Zj (6) 

where j G T- The ag^r, agj, and a^j are channel gains 
between the nodes (see Fig. |^. We later relate these gains to 
distances between the nodes. The Zr and Zj, j = 1,2,..., N, 
are independent and identically distributed Gaussian random 
variables with zero mean and unit variance. We may alterna¬ 
tively write © and (|^ in vector form as 

Y, = A,X -k Z, (7) 

where X = {Xg XrY, Yj = {Yr Yj^, Z = {Zr ZjY, and 



III. Quasi-Concavity in SNRs and p 
A. Cut-Set Bound 

We consider two characterizations of Res- First, let aj = 
(osj Crj) be the second row of Aj, let Qx be the covariance 
matrix of X (see Appendix and let det M be the deter¬ 
minant of the square matrix M. The cut-set rate ( [T^ can be 
expressed as the maximum of 


-Rcs(Qx) 


mm mm 

^<3<n 


log (aJ Qx -k l) , 


Y / det Q^y7 


XrY 


( 16 ) 


over the convex set of Qx with diagonal entries Pg and Pr- 
The first logarithm in ( fTfi] ) is clearly concave in Qx. The 
second logarithm is concave in Q(yt’ Xr-Y Appendix |B|) 
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and Q(yT x^)t is linear in Qx- To prove the latter claim, 
observe tfiat 

Q(YJ Xr-)^ = AjQxAJ ^ q ^ 

where Aj = [0 1]^^ and I 2 is the 2x2 identity 

matrix. Hence i?cs(Qx) is concave in (the convex set of) 
Qx because it is the minimum of 2N concave functions. 

Suppose next that we wish to consider p and the SNRs 
individually rather than via Qx. Dehne the vector 


B. DF Rate 

Consider the functions 

5*(p,S) = (l-p2)SNR,., (23) 

Rdf{p,S) = ^min^min (C(/j(p,S)),C(p*(p,S))) . (24) 

Theorem 2: S) is concave in p, concave in S, and 

quasi-concave in (p^,S) for 0 < < 1 and non-negative S. 

Proof: The proof is similar to that of Theorem H ■ 
Corollary 2: Rdp{p,S{P)) is quasi-concave in (;(r, P) for 
0 < < 1 and non-negative P. 

Proof: See the proof of Corollary ■ 


S = (SNR,,„SNR,,i,-.. ,SNR«,jv,SNR,,i,--- ,SNR,,jv) 
and the functions 


C. DT Rate 

The DT rate ([T3|l is clearly concave in S and P. 


fjip,S) = SNRsj + SNR^j -f 2p^SNRsjSNR^j (18) 
p, (p, S) = (1 - p2) (SNR,,, + SNR,,,) (19) 

Rcs{p,S) = min min (C(fj(p,S)),C(gj(p,S))). (20) 


We establish the following results. 

Lemma 1: fj{p,S) and p, (p, S) are concave in p, concave 
in S, and quasi-concave in (p^, S) for 0 < p^ < 1 and non¬ 
negative S. 

Proof: Concavity with respect to p is established by 
observing that /, (p, S) is linear in p, and pj(p, S) is linear 
in —p^ which is concave in p. 

Consider next concavity with respect to S. The Hessian of 
fj{p,S) with respect to S has only one non-zero eigenvalue 


P 

2 


SNR2 , -fSNR^, 
SNR^JSNR^// ■ 


( 21 ) 


D. QF Rate 


Consider the functions 


h,(S) = SNR,,, + 


i?QF(S) 


min 


SNR^ jSNRg 


SNRsj + SNR.pj + SNRs,r + 1 

C(h,(S)). 


(25) 

(26) 


We establish the following results. 

Lemma 2: hj{S) is quasi-concave in (SNR,,, SNR,,,) for 
non-negative SNRs. 

Proof: Substitute a = SNR,,,-, b = SNR,,,, k = SNR,,, + 
1 into the second function of Lemma and apply Lemma 
Part [2 _ ■ 

Theorem 3: Rqp{S) is quasi-concave in the non-negative 
SNRs if the SNR,,,, j = 1, 2,..., n, are held hxed. 

Proof: Apply Lemmas and Lemma Parts and ■ 


Thus, fj{p, S) is concave in S. The function p,(p, S) is linear 
in S, and thus concave in S. 

Now consider quasi-concavity with respect to (p^, S). Sub¬ 
stituting a — SNR^, 6 = SNR,,,c = p^ into the hfth 
function of Lemma we hnd that fj{p,S) is quasi-concave 
in (p^, S). For the p, (p, S), observe that ab is quasi-concave 
for non-negative (a, b) (see the hrst function of Lemma |^. 
This implies 

(Aoi -I- Aa2)(A6i -I- A&2) > min (ai 5 i, 0262) ( 22 ) 

for 0 < A < 1, and where A = 1 — A. Substituting Oi = 1 — pj 
and bi = SNR,,,,.i-|-SNRs,,,i for j = 1, 2, we hnd that p, (p, S) 
is quasi-concave in (p^,S). ■ 

Theorem 1: Rcs{p, S) is concave in p, concave in S, and 
quasi-concave in (p^, S) for 0 < p^ < 1 and non-negative S. 

Proof: Rcsip, S) involves taking logarithms and minima 
of (quasi-) concave functions. The results thus follow by 
applying Lemma [2 and Lemma Parts and |2 ■ 

Corollary 1: Consider S as a function of the powers P = 
{Ps,Pr)- Then i?c's(p, S(P)) is quasi-concave in (p^,P) for 
0 < p^ < 1 and non-negative P. 

Proof: The proof follows from Theorem [2 and because 
S is a linear function of P. ■ 


IV. Quasi-Concavity in Relay Position 
Suppose the channel gain for the node pair {i,j) is 


(lip 



(27) 


where , is a “fading” gain, , = ||i —j|| is the Euclidean 
distance between the positions i and j of nodes i and j, 
respectively, and a > 2 is a path-loss exponent. We thus have 


SNR,,, 





We establish quasi-concavity results in p^ and r. 


A. Cut-Set Bound 

Consider the functions ([T§-@ but relabeled as /,(p, r), 
p,(p, r), and Rcs{p,t:) to emphasize the dependence on the 
considered parameters. 

Lemma 3: /,(p, r) and gj{p,r) are quasi-concave in r for 
hxed p, 0<p< 1. Furthermore, /, (p, r) is quasi-concave in 
(p^, r) for 0 < p^ < 1. 

Proof: Consider the functions 


fMDn = 


^sjPs ^r,jPr 




D° 


2p, 


' ^spPs ^rjPr 


D° 




(28) 
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= + ( 29 ) 

which are quasi-linear in for fixed p since they are 
decreasing in D“. But ^ and are convex functions 
of r for a > 1, and thus Lemma Part establishes that 
fj {p, r) and pj {p, r) are quasi-concave in r for fixed p. 

Next, substitute a = and b = ^ into the third function 
of Lemma|^ and use Lemma|^ Part|T[ to show that fj {p, _D“) 
is quasi-concave in But fj is decreasing in and 

Dfj is convex in r, so Lemma Part establishes that 
fj{p,r) is quasi-concave in {p^,ry ■ 

Unfortunately, pj is quasi-convex (and not quasi-concave) 
in {p^,D°‘). To see this, substitute a = and b = p^ into 
the fourth function of Lemma Quasi-concavity would have 
been useful since it would have permitted using Lemma 
Parts 1^ and to establish the quasi-concavity of 


-Rcs(i’) = max 

p 


min min 

l<i<7V 




(30) 


Based on our numerical results, i?c's(r) seems to be quasi¬ 
concave in r but we have been unable to prove this. Never¬ 
theless, Lemma suffices to establish an intermediate result 
which is useful in Sec. [V] when we study p = 0. 

Theorem 4: Res {p, r) is quasi-concave in r for fixed p, 

0 < p < 1. 

Proof: Rcs{Pj^) is the minimum of functions that are 
quasi-concave in r. Lemma Part thus establishes the 
theorem. ■ 


B. DF Rate 

The quasi-convexity of pj{p,D°‘) relaxes for the DF 
rate Consider the negative of the fourth function of 

Lemma 1^ with fci = 0: 

f{a,b) = {l-b)k 2 /a. (31) 


This function is quasi-linear in (a, b) since both its superlevel 
and sublevel sets are convex. This result implies the follow¬ 
ing theorem. We again consider the functions (|2^-(|24ll but 
relabeled as p*{p,v) and RDp{p,r). We further define 


pNp,D“) = (l-p^) 


2 \ ^s,rPs 




(32) 


i?DF(r) = max 

p 


imm^min {C{fj{p, r)), C{p*{p, r))) 


(33) 


Theorem 5: Rpp{p,r) is quasi-concave in (p^,r) for 0 < 
p^ < 1, and Rp)p{r) is quasi-concave in r. 

Proof: Pj{p,D°‘) is quasi-linear in and de¬ 

creasing in Furthermore, is convex in r, and thus 
Lemma Partj^ establishes that p*(p, r) is quasi-concave in 
(p^,r). Rp,p{p,r) is therefore quasi-concave in r, as it is the 
minimum of quasi-concave functions (see Lemma ^ Part |^. 
Furthermore, Rp,p{r) is concave in r by Lemma |8 Part|^ ■ 



V. DF Performance 

This section presents numerical results for the DF strategy 
and compares them to results from our previous work 0-0. 
We consider 1-, 2-, and 3-dimensional MRCs with different 
numbers N of destination nodes. For simplicity, we consider 
the low-SNR or broadband regime where 

C(SNR) = ^log(l + SNR) ^ isNR. (34) 

In other words, we consider the cut-set and DF rates without 
the logarithms. This approach is valid not only in the limit of 
low SNR, but more generally because we proved our quasi¬ 
concavity results without taking logarithms. Furthermore, in 
the low-SNR regime the rates of full-duplex and half-duplex 
transmission are the same under a block power constraint. 

We choose Pg = Pr = P = 1, a = 2, and = 1 for all 
node pairs {u,v). We study both coherent transmission where 
p is optimized and non-coherent transmission with p = 0. 
The rates are in nats/channel use. Alternatively, suppose we 
use sine pulses sampled at 2W samples per second, where W 
is the (one-sided) signal bandwidth. Suppose further that the 
(one-sided) noise power spectral density is 1 Watt/Hz. Then 
at low SNR the rates in nats/channel use are the same as the 
rates in nats/sec. 


A. One Dimension 


Consider a relay channel {N = 1) where the source is at 
the origin (s = 0) and the destination is at point 1 (1 = 1). 
Fig. shows the low-SNR cut-set bounds, DF rates, and the 
routing-based DF (RDF) rates developed in ||^. For example, 
the non-coherent rates with p = 0 are 


R 


DF 


Is- 111 


P 

Rrdf max — 
0</9<l 2 


, ^s,r 

||r- 1|1“’ ||s-r 

Cr,l fs, 

Ir- 1 


l|s-l||“ 


(35) 


( 36 ) 
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Observe that all curves are quasi-concave (but not concave) in 
r. Theorems 1^ and [^predict the quasi-concavity for all curves 
except for the coherent cut-set bound. 

The best relay positions for the two strategies are different. 
For example, r = 0.5 maximizes Rrdf while the r max¬ 
imizing Rdf is closer to the source. This is because when 
the source transmits, the relay and the destination listen and 
the destination “collects” information. The relay can thus be 
positioned closer to the source while maintaining the same 
information rate from the source to the relay, and from the 
source-relay pair to the destination. At the optimal positions, 
we compute Rdf ~ 2.25P nats/sec and Rrdf — 2P 
nats/sec, so the DF gain is « 12%. We remark that the non¬ 
coherent cut-set bound coincides with (the non-coherent) Rdf 
and Rrdf for a large range of r. 


B. Two Dimensions 


Consider N = 5 destinations positioned on a square in the 
two-dimensional Euclidean plane with the source node at the 


origin. Fig. 4a plots the non-coherent Rdf as a function of 
the relay position. The best relay position is shown by a 
circle and the corresponding rate is Rdf ~ 0.01 IP nats/sec. 
Fig. 1^ plots the low-SNR two-hop rate 

^s,r ^r,j 


P 

R 2 H mm — mm 
l<j <5 2 


s - r 




(37) 


as a function of the relay position. The best relay position 
r 2 D is shown by a circle and the corresponding two-hop rate 
is R 2 H ~ 0.0095P nats/sec. The non-coherent DF gain is thus 
« 12 %. 


Fig. 4b and Fig. 4d show contour plots for Rdf and R2h- 


The contours form convex regions, as predicted by Theorem]^ 
Again, the relay position maximizing Rdf lies closer to the 
source than the relay position maximizing R 2 h- 


VII. Conclusion 

Various quasi-concavity results were established for AWGN 
MRCs. Quasi-concavity lets one use efficient optimization 
algorithms to maximize the rate and to place the relay. 

Appendix A 

Concave and Quasi-Concave Functions 
A. Definitions 

Consider the following sets. The domain Df of a real-valued 
function / : R" —>• R is the set of arguments for which / is 
dehned. The hypograph and hypergraph of / are the respective 

"W/= {(x,y) I 2 /</(x)} (38) 

"W/= {(x,y) I 2 />/(x)}. (39) 

The superlevel and sublevel sets of / with respect to /3 S R 
are the respective 

Sf,fi = {x|/(x) > /?} (40) 

= {x|/(x) < /?}. (41) 

Concave and quasi-concave functions can be dehned via the 
convexity of these sets. Recall that a set S, S C R", is convex 
if for any two points xi and X 2 in S and for any A satisfying 
0 < A < 1 we have 

Axi -f Ax2 S S (42) 

where A = 1 — A. Suppose that Df is convex. The function 
/ is concave over Df if and only if its hypograph Hf is 
convex. Similarly, / is convex over Df if and only if Rf is 
convex. The function / is quasi-concave over Df if and only 
if all its superlevel sets are convex, and / is quasi-convex over 
Df if and only if all its sublevel sets are convex. A function 
that is quasi-convex and quasi-concave is called quasi-linear. 
For example, any non-increasing or non-decreasing function 
is quasi-linear. 


C. Three Dimensions 

Consider N = 5 destinations positioned in 3-dimensional 
Euclidean space as in Eig. The hgure also shows the convex 
hull (a polyhedron) of the points. The points r^p and r^jj 
denote the relay positions that maximize the non-coherent 
Rdf and R2H, respectively. We remark that rj^p and r^jj 
remain unchanged if more destinations are positioned inside 
the polyhedron. This is because the points in the polyhedron 
can receive at least the same rate as the worst of the hve nodes 
at the corner points. 

VI. Discussion 

We considered real AWGN channels. Eor complex AWGN 
channels, quasi-concavity in r will likely not be valid because 
the phases of the channel gains a„ „ will change with r, and 
this means that p must be adjusted differently for different 
destination nodes. We remark that this effect is “local” in the 
sense that for large carrier frequencies the phase variations are 
sensitive to changes in r. A pragmatic approach would then 
be to optimize r for non-coherent transmission (p = 0) even 
if beam-forming is permitted. 


B. Basic Properties 

Two properties of concave and quasi-concave functions are 
as follows; these properties are often used as the dehnitions of 
such functions. Similar properties exist for convex and quasi- 
convex functions. 

Lemma 4: The function / is concave if and only if 

/(Axi -I- Axa) > A/(xi) -f A/(x 2 ) (43) 

for all Xi and X 2 in I?/ and for all 0 < A < 1. 

Lemma 5: The function / is quasi-concave if and only if 

/(Axi + AX 2 ) > min(/(xi), /(X 2 )) (44) 

for all Xi and X 2 in D/ and for all 0 < A < 1. 

The next two properties assume that / is twice differentiable 
and that Df is convex. Let H f (x) and B/ (x) be the respective 
Hessian and bordered Hessian of / at x. 

Lemma 6: (see ||^ Sec. 3.1.4]) / is concave if and only if 
H/ (x) is negative semidehnite for all x G Dy. 

Lemma 7: (see p. 771]) / is quasi-concave on the open 
and convex set Df if the determinants D2 , , • ■ •, Dn of the 

respective second to nth leading principal minors of Bj(x) 
satisfy < 0 for fc = 2, 3,..., n and for all x G Df. 
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Contours of the function Rbf 


X - axis (meters) 

(b) 



■ 10-2 


y (meters) x (meters) 

(C) 

Fig. 4; (a) Rdf for N = 5; (b) Contour plot of Rdf', (c) R 2 H 


X - axis (meters) 

(d) 

for the same network; (d) Contour 


plot of R 2 H- 


Contours of the function R 2 H 


2 



Fig. 5: N = 5 destination geometry in three dimensions. 
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C. Compositions Preserving Quasi-Concavity 


4) We have D 2 < 0 and D 3 > 0 for 


The following compositions preserve quasi-concavity. 
Lemma 8: Suppose / and fi, 1 < j < n, are quasi-concave, 
then so are the functions 

1 ) h = kif -f k 2 , where ki > 0 and k 2 S R; 

2 ) h = min f, ; 

3) h = g o f where / is quasi-concave and g is non¬ 
decreasing; 

4) h{a) = supbgg /(a, b) where S is a convex set; 

5) h{a,h) = f{g{a),h) where g is convex and /(a, b) is 
non-increasing in a for fixed b. 

Proof: Properties l)-4) are standard (see |[^ Sec. 3.4]). 
For property 5), observe that 

h{Xai -\- Aa 2 , Abi -I- Ab 2 ) 

= + ^^ 2 )! ^bi -f Ab2) 

(a) 

> /(A5(ai) -f Ap(a2), Abi -f Ab2) 

(b) 

> min(/( 5 (ai),bi),/(g(a 2 ),b 2 )) (45) 

where (a) follows because g{Xai -f Aa 2 ) < Ap(ai) -f Ap(a 2 ) 
and /(a, b) is non-increasing in a. Step (b) follows because 
/ is quasi-concave. ■ 


D. Examples of Quasi-Concave Functions 


We establish quasi-concavity for several useful functions. 

Lemma 9: The following functions are quasi-concave for 
X = (a b c) with non-negative entries. 

1 ) /(x) = ab 

2 ) /(x) = for a positive constant k 

3) /(x) = ki/a -\- 2 y/k 2 b/a for positive constants ki, k 2 

4) /(x) = —{l — h){ki-\-k 2 /a) for positive constants ki,k 2 

5) /(x) =0 4-6-1- 2s/abc 

Proof: We use bordered Hessians B/ (x) and the deriva¬ 
tives Dk of their fcth leading principal minors, k = 2,3,... ,n. 
The results extend to non-negative x by using continuity at 
zero values, except for the third and fourth functions where 
0 = 0 makes the functions undefined. 


1) We have D 2 < 0 and D 3 > 0 for 


B/(x) = 


0 be ac 
be 0 c 
oc c 0 


2) We have D 2 < 0 and > 0 for 


B/(x) = 


0 

b{b-\-k) 

{a-\-b-\-k)^ 


b{b-\-k) 

(a+b+A:)2 

-2b(b-\-k) 

(a-\-b-\-k)^ 


a{a-\-k) 2ab-\-{a-\-b-\-k)k 

\ {a-\-b-\-k)‘^ (a+6+fc)3 


a{a-\-k) 

(a+fc+fc )2 

2ab-\-{a-\-b-\-k)k 

{a-\-b-\-k)^ 

— 2a{a-\-k) 


3) We have D 2 < 0 and D 3 > 0 for 


B/(x) = 


0 


fci + 'N/fc2Q^ 


ki-\--\/k2ab 4fci+3\/A:2a^ 

2 a^ 


2a3/2^ 


\ 

ab 

•/kg 

263/2 


5 ) 


B/(x) = 


0 

(l-6)fc2 

a 2 

ka 


(l-b)k2 

a3 

2(l-6)fe2 

a3 

_ hi. 


fci + f 

k 2 


+ ^ 0 

We have D 2 < 0, > 0 and £>4 < 0 for 


\ 

/ 



Appendix B 

Covariance Matrices and Concavity 

The covariance matrix of a real-valued random column 
vector V is 


Qv = E[(V-E[V])(V-E[V])^] . (47) 

A useful property of covariance matrices is as follows (see m 
p. 684]). If is a principal minor of Qv, then the function 


/(Qv) = log 


det Qv 
det 


(48) 


is concave in Qv- 
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